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Abstract 

Baker (2008) introduced a new class of bivariate distributions based 
on distributions of order statistics from two independent samples of size 
n. Lin-Huang (2010) discovered an important property of Baker's dis- 
tribution and showed that the Pearson's correlation coefficient for this 
distribution converges to maximum attainable value, i.e. the correlation 
coefficient of the Frechet upper bound, as n increases to infinity. Bairamov 
and Bayramoglu (2011) investigated a new class of bivariate distributions 
constructed by using Baker's model and distributions of order statistics 
from dependent random variables, allowing high correlation than that of 
Baker's distribution. In this paper a new class of Baker's type bivariate 
distributions with high correlation are constructed on the base of distri- 
butions of order statistics by using an arbitrary continuous copula instead 
of the product copula. 

Keywords: Bivariate distribution function, FGM distributions, cop- 
ula, positive quadrant dependent, negative quadrant dependent, order 
statistics, Pearson's correlation coefficient. 



1 Introduction 



Huang and Kotz (1999) introduced new modifications of classical Farlie-Gumbel- 
Morgenstern (FGM) distribution introducing additional parameters. The new 
Huang-Kotz FGM distributions allow correlation fiigfier than the classical FGM 
and because of the simple analytical form aroused interest of many researchers. 
Last years there appeared many papers dealing with the modifications of FGM 
distribution allowing high correlation. For a related works on this subject see 
Bairamov et al. (2001), Amblard and Girard (2002), Bairamov and Kotz (2002), 
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(2003), Fisher and Klein (2007) among others. Baker (1998) use a novel ap- 
proach connected with the FGM distribution and introduce a new class of bivari- 
ate distributions based on the distributions of order statistics. Baker (2008) con- 
siders independent random variables Xi, Yi, from two univariate distributions 
with distribution functions (c.d.f.) Fx and Fy, respectively. The correspond- 
ing probability density functions (p.d.f.) are fx and /y. Let Ui = m\n(Xi,X2) 
and U2 = max(Xi, X2), Vi = min(yi, F2) and V2 = max(Yi, F2)- To obtain pos- 
itive correlation Baker randomly chooses either the pair Ui,Vi or (72, V2. The 
random numbers are now positively correlated, but the marginal distributions 
are unchanged, because the random choice of either of two order statistics from 
a distribution gives a random variable from that distribution. To obtain a neg- 
ative correlation, either Ui,V2 or U2,Vi chosen. The bivariate distribution of a 
randomly chosen pair of order statistics is 

1/2 X {Fr{x)Fr{y) + F}i\x)Fl.-\y)} 

and on choosing either a pair of order statistics with probability q or the 
original independent random variables X, Y with probability 1 — g, the bivariate 
distribution function is 

H{x,y) = {l-q)Fx{x)Fy{y) + {q/2){Fr{x)Fr{y) + F}i^{x)F^-\y)] 
= Fx{x)FY{y){l + q{l - Fx{x)){l - Fyiy)}. 

In general let Xi , X2, X„ and Yi , 1^2, ■■■Yn are independent and identically dis- 
tributed (i.i.d.) random variables with distribution functions (d.f.) Fx and Fy , 
respectively. Let Xfc.„ and Y^^.^, k = l,2,...,ri be corresponding order statis- 
tics and i^f "(x) = P{Xk:n < x}, Fp'^iy) = P{Yk:n < y}- Baker's bivariate 
distribution function is now defined as 

n 
k=l 

H[^\x,y) = -J2F^--{x)F^-'^+'---{y). (2) 

k=l 

For Baker's bivariate distribution iJ^""-* (a;, y) with exponential marginals Fx (x) = 
Fy{x) = 1 — e~^,x > 0, the Pearson's correlation coefficient is p„ = 1 — 
n Sfc=i i' which increases monotonely to 1. 

As a generalization of ([1]) and ([5]) Baker also introduces 

n n 

i7(")(x,2/) =EEr,,F|"(x)F^"(y), (3) 
fc=i 1=1 

where rki > and J^k^i'^ki = Z^Li ^fc' = h ^^"^ ^ l,2,...,n. Lin 

and Huang (2010) proved that ([3]) does not contain members with a corre- 
lation higher than that of ([1]). That is why the best bivariate distribution 
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with higher positive correlation among the members ([3]) is (H)). Similar con- 
sideration holds true for the negative correlation. Lin and Huang (2010) in- 
vestigate the conditions under which the correlation for ([T]) converges to the 
limit. In particular, they show that if either (i) X > b,Y > c a.s. for some 
6,c e M and E{Xk:n) > F^\^) and E{Yu.,r,) > Fy'i^) for aU (fc,n), 
or (ii) X < b,Y < c a.s. for some 6, c G M and E{Xk:n) < Fx^in^ ^'^'^ 
E{Yk;n) < Fy^{—) for all {k,n) then lim p„ = p*, where p* is the correlation 

coefficient of the Frechet-Hoeffding upper bound H+{x, y) = min(i^x{a;), Fyiy)) 
(see Frecheet (1940)). The results presented in the paper of Lin and Huang 
(2010) makes Baker's distribution attractive. 

Recently, Bairamov and Bayramoglu (2011) observed that if in the Baker's 
model instead of independent random variables one uses the dependent random 
variables (X, Y) with positive quadrant dependent (PQD) joint distribution 
function F{x,y), then the correlation increases and for negative quadrant de- 
pendent F{x, y), decreases. More precisely, let {Xi,Yi), {X2, 12), (X„, 1^) be 
a bivariate sample with joint distribution function F{x,y) = C {Fx (x) , Fy (y)) . 
Bairamov and Bayramoglu (2011) consider the following bivariate distribution 
functions constructed on the basis of the Baker's idea: 



1 

K^+\x,y) = -y^P{Xr:n<X,Yr..n<y} 
r=l 

K^"\x,y) = -y^P{Xr:n<X,Yn-r+l:n<y}, 

71 ' * 



n 

r=l 



where Xi-n and Yj:„ are the ith and jth order statistics constructed on the 
basis of bivariate observations {Xi,Yi), {i — l,2,...,n) with joint distribution 
function F{x^ y) = P{Xi < x,Yi < y} and marginal distribution functions 
Fxix) = F{x, 00), Fviy) = F(oo, y) so that Xi:„ < X2:n < ■ ■ < X^-.n, < 
5^2:n < • • • < Yn:n- The joint d.f. of Xr:n and Ys:n IS given in David (1981) (see 
also Arnold et al. (1992)) as 

P{Xr:n < X, Y^-.n < y} 

71 n b 

i—r j—s k—a 

where 

n' 

^("'fe'''^')= fc!(,_fc),(j_fc),(„„,_j + fc)! ' « = niax(0,z-fj-n), 6 = min(^,J) 

(5) 



and 



Pll = F{x,y) 

P12 = Fx{x)-F{x,y) 

P21 = FY{y)-F{x,y) 

P22 = l~Fx{x)-FY{y)+F{x,y). (6) 
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Then 



^ n 71 n b 

r— 1 i—r j—r k—a 
^ n 71 71 b 

r— 1 z— r J— n— r+1 k—a 

It is clear that the marginal distributions of K'^\x,y) and K^\x,y) are 

again Fx(a^) and Fy^y), i.e. oo) — Fx{x) and i^i"^(oo, y) = Fyiy). It is 

shown that for a PQD joint distribution function F{x, y) the positive correlation 
of K^\x,y) is higher than that of H^\x,y) and for NQD F{x,y) and the 

negative correlation of K^\x,y) is smaller than that of H'^^\x,y). 

In this note we consider a new class of bivariate distribution functions using 
Baker's construction and considering any copula C{u,v) instead of product 
copula C(w, v) = n(M, v) = uv. It follows that for this new class of distributions 
if C{u,v) is PQD, i.e. C{u,v) > uv, for all {u,v) £ [0,1]^, then the Pearson's 
correlation coefficient is higher than that of Baker's distribution. Similarly, if 
the copula is NQD, i.e. C{u,v) < uv, for all {u,v) G [0,1]^ then then the 
Pearson's correlation coefficient is smaller than that of Baker's distribution. All 
modifications constructed on the Baker's idea we call as Baker-Lin-Huang Type 
distributions. The distribution ([7]), (H)) we call Baker's Type I BB and the new 
distributions introduced in this paper Baker's Type II BB distributions. 

2 New bivariate Baker's Type II BB distribu- 
tions based on a copula approach 

Let Xi, X2, Xn and Yi,Y2, ...Yn be i.i.d. random variables with d.f.'s Fx 
and Fy, respectively. Let Xfe.„ and Yfc:„, k — 1,2,..., n be corresponding order 
statistics and F^-"{x) = P{Xk..n < x}, F^--"{y) = P{Yk..n < y}- 

Recall that a two-dimensional copula is a function C(x,y) from [0,1]^ = 
[0, 1] X [0, 1] to [0, 1] with the properties: 

1. C{x, 0) = = C(0, y), C{x, l)=x and C(l, y) = y; 

2. For every xi,X2, 2/1,2/2 such that < a:i < X2 < 1 and < ?/i < 2/2 < 1 

C{x2,y2) - C{x2,yi) - C{xi,y2) + C{xi,yi) > 0. 

According to Sklar's Theorem if Fx^vix, y) is a joint distribution function with 
continuous marginal distributions Fx{x) and Fy{y), then there exists a unique 
copula C such that F{x,y) = C {Fx (x) , Fy (y)) . Theory and applications of cop- 
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ulas are well documented in Nelsen (2005). Let C{u, v) be any copula. Consider 

1 " 

G^^\x,y) = -Y,C{F^^-{x),F^--{y)) (9) 

GL"^(x,y) = -Y,C{F^^-{x),F-~'^+'---{y)). (10) 

^ k=l 

From the properties of a copula it follows that the marginal distributions of 
G^^\x,y) and G^^\x,y) are Fx and Fy, respectively. In fact, 

gV"^(^,^) - -Y^C(F^x''{x)^) = -Y^F';^-{x)^Fx{x) 

k=l k=l 

-.71 1 

GV"^(^,y) - -Y,C{l,F^-{y))^-Y,F^--{y)^Fy{y). 

k=l k=l 

Similarly, g'"^\x,oo) — Fx{x) and G^I^\oo,y) = Fyiy). Hereafter, we will de- 
note by the Pearson's correlation coefficient between any random variables 
X and Y with joint distribution function H(x, y). It is clear that if G{u, v) = 
n(u,w) = uv then G^]^\x,y) — H^^\x,y) and G^^\x,y) — H^^^\x,y). Since 
C {Fx" (x) , Fy" (y)) is a bivariate c.d.f. (with marginals F^'^ix) and Fy^iy)) 
then G^^\x,y) is obviously a bivariate c.d.f. as a convex combination of bi- 
variate c.d.f. 's. The copula used in construction © and ([TU]) will be called the 
"kernel" copula for Baker's Type II BB distribution. 

Theorem 1. If C{u,v) is PQD then p„(n) > p„(n) and if G{u,v) is NQD 
then pg(„) < 

Proof. Since C(u,i;) > to for aU (u, ^;) G [0, 1]^, then C(Ff "(x), F^^"(y)) > 
F^-"{x)FY"{y) for all {x,y) e M^. From the Hoeffding's formula (see Hoeffding 
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(1940)) for correlation coefficient one has 



oo oo 



= —=!===[ [ [d^\x,y)^Fx{x)FY{y)]dxdy 
+ \JVar{x)V ar{Y) ' ' 



1 



-OO — OO 

oo oo 



n 



— oo — oo 



n 
fc=l 



^Var{x)Var{Y) 
-Fx{x)FY{y)]dxdy (11) 



1 



yJVar{x)Var{Y) 
1 



[H^^\x,y)-Fx{x)FY{y)]dxdy 



— OO — OO 

oo oo 



^Var{x)Var{Y) 
-Fxix)FY{y)]dxdy (12) 
and > . Similarly, Pg,(„) < . 



+ + 



Example 1. Baker's Type II BB distributions with FGM "kernel" 
copula. Let 

C{u,v) = Mw(l + a(l -u)(l - v)),{u,v) e [0, < a < 1. (13) 
Consider 

G^^\x,y) = lf]Fl-(x)F^-(y)(l+ (14) 
ail-F^"ix))il-F^-"iy)), < a<l. 

n 

G^^\x,y) = -Y,F'^-{x)F--^+'---{y){l + 

k=l 

a{l-F^-"{x)){l^Fr''+'--'^{y)), -1 < «<0. (15) 



Let Fx{x) = x,0 < X < 1 and (y) = y, < y < 1. Since the FGM 
copula (IT^ is PQD for a > and is NQD for a < 0, then > p„(") for 

> 0, /OqC") < /OjLf("), for a < 0. It is clear that if a = then G^^\x,y) = 

H^\x,y) and G^J^\x,y) — H^_^\x,y). In the following table we present some 
numerical values of p^(n) , p^cn) , , p^cn) , P^C") and Pjj(n) for Uniform{0, 1) 

marginals. The numerical calculations are made in MATLAB which is one 
of the commonly accepted packages for coding mathematical models since it 
has built-in functions for probability distributions and allows probabilistic and 
mathematical operations. 
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n 


2 


4 


6 


8 


10 


12 


15 


20 


Pa'f 


0.5467 


0.7258 


0.8039 


0.8475 


0.8753 


0.8945 


0.9144 


0.9348 




0.5133 


0.6915 


0.7761 


0.8247 


0.8561 


0.8779 


0.9006 


0.9241 




0.3333 


0.6000 


0.7143 


0.7778 


0.8182 


0.8462 


0.8750 


0.9048 


PaM 


-0.5467 


-0.7258 


-0.8039 


-0.8475 


-0.8753 


-0.8945 


-0.9144 


-0.9348 


PkM 


-0.5133 


-0.6915 


-0.7761 


-0.8247 


-0.8561 


-0.8779 


-0.9006 


-0.9241 


Ph^"-) 


-0.3333 


-0.6000 


-0.7143 


-0.7778 


-0.8182 


-0.8462 


-0.8750 


-0.9048 



Table 1. Correlation coefficients and with 

F{x,y) = xy{l + a{l - x){l - y)), a ^ l,Fx{x) = x,FY{y) = y,0 < x,y <1 
and p^(^) , Pj^(^) and p^{r^) , p^^m with Uniform{0, 1) marginals. 



It can be observed from the Table 1 that p^M < Pj^(n) < Pjj(n) < PjjM < 



In the following table the values of the correlation coefficients for p„(n) , yO^(»i)/9„(n) , 
Pg(„),/9^(n) and pjjM for Uniform{0,l) and Exponential{l) marginal distri- 
butions. 



n 


2 


4 


6 


8 


10 


12 


15 


20 


Pc^f 


0.4811 


0.6343 


0.7001 


0.7367 


0.7600 


0.7762 


0.7929 


0.8102 


PK'f 


0.4426 


0.5951 


0.6682 


0.7105 


0.7379 


0.7572 


0.7772 


0.7980 


P^i^) 


0.2886 


0.5196 


0.6185 


0.6735 


0.7085 


0.7327 


0.7577 


0.7835 




-0.4811 


-0.6343 


-0.7001 


-0.7367 


-0.7600 


-0.7762 


-0.7929 


-0.8102 




-0.4426 


-0.5951 


-0.6682 


-0.7105 


-0.7379 


-0.7572 


-0.7772 


-0.7980 




-0.2886 


-0.5196 


-0.6185 


-0.6735 


-0.7085 


-0.7327 


-0.7577 


-0.7835 



Table 2. Correlation coefficients pQ{n) and p^in) with 

F{x,y) = xy{l + a{l - x)(l - y)),a = l,Fx{x) = x,0 < .t < l,Fy(y) = 1 - exp(-y),2/ > 
and Prf(n), Prf(n) 0,11(1 wlth thc Same marginals. 



Again, from the Table 2 we have p^^n) < Pj^(n) < Pj^M < PjjM < Pj^M < 
Pain)- 
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3 Copula representation of Baker's Type II BB 
distribution 

Denote by n(i, s) — ts, {t, s) G [0, 1]^, a product copula. Let {Xi,Yi), i ~ 1,2, ..n 
be a bivariate sample with joint distribution function F(x, y) = C{Fx{x), Fyiy))- 
Consider the joint distribution of order statistics Xr.n and Yg-.n given in Q: 

Fx^:^.Y^.Ax,y) = P{Xr:n < X,Y,.,n < y} 

n 71 b 

= EEE^("'fc;*'-?')^'n^^i2V2r'pr"''"', (16) 

i—r j—s k—a 

where c{n,k;i,j), pii,pi2,P2i,P22 are given in (O,®. The copula of vivari- 
ate distribution Fx^.^^Ysni^iV) presents an interest. Denote this copula as 
Cr,s:n{t, s) , then 

Fx^:,.yr-.Sx.y) = Cr.s:n{Fx^.Ax),FY,Ay)). (17) 
It is well known (David (1981)) 

Fx..Ax) = J2(])pxi^^(^-Fxix)r-^ 

Fx{x) 

U-'\l - uT-^du - l„r-r+l{Fx{x)) 



B{r, n — r + 1) 





and 



FY..Ay) - E( J^^(2^)(i-^^(2/))^ 

FY{y) 



^ • u'-\l-ur-'du = Is,n-s+l{FY{y)), 



B{s,n- s + 1) 



p 

where Ia.b{p) = B(a,b) j ^"""'^(1 — u)'^~^du is an incomplete Beta function. De- 



note by I^l{p) the inverse of Ia.b{p)- Let Fx^.^{x) — Ir.n-r+iiFx{x)) — t and 
Fy^.Jjj) = /.,„-.+i(Fy(y)) = .s. Then 

x = F^\l-;^_^^,iFxm and y = Fy\l-^^_,^,{Fy{s))). (18) 

From ini) and ^ one has 

a,s:„(t,s) = Fx_,y_ {Fx\lr,Lr+lit)),Fy\l-l_,+ ,{s))) . (19) 
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Therefore a copula of joint distribution of order statistics Xr-.n and Ys-.n is 
Cr,s:n(t, s) givcu in (|19p . In a special case if r = s = n one has 



(20) 



since /„,i(i) = and /,7_l(w) = w^/". 

Analogously, one obtains the copula of joint distribution of order statistics 
Xi:n and Yi-n as 



+^^(F^i(l-(l-t)i/"),Fy(l-(l-.)i 
by noting that the joint distribution function of Xi.„ and Yi;n is 



(21) 



Fxi.,„,Yt,„{x,y) = P{Xi.,n < x,Yx;n < y} 

= (1 - (1 - Fx{x)r) + (!-(!- FY{y)r) - 1 + F"(x,y) 
= (1 - (1 - Fx{x)r) + (!-(!- Fy(y))") - 1 

+ {1 - Fxix) - Friy) + Fix,y)r 
= Ci.n:n{l-{l-Fx{x))r,l-(l-FY{y)r)- (22) 

3.1 The Baker's Type II BB distributions with "kernel" 
copula of bivariate FGM extreme order statistics 

If the underlying distribution is classical FGM, i.e. 

Fix, y) = Fx{x)FY{y){l + a{l - Fx{x)){l - Fyiy))), -!<«<!. (23) 
one obtains from (PU)) 

-in 

C„:„(i, s) ^ts\l + a{l - <i/")(l - si/")J (24) 

and from (HH) 

Cl:„(i,s) = t + S-l+ [(l-i)'/" + (l -«)'/"-! 

+ (l-(l-t)i/")(l-(l-s)i/") 

x{l + a(l-t)i/"(l-s)i/"} . (25) 



It follows that lini C„:„(f, s) — ts — Il{t,s) and lim Ci:„(i, s) = n(i, s). 
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since lim [l + a(l - - s^/")] " = 1 and 



lim (1 - t)i/" + (1 - s)i/" - 1 

+ (1 - (1 - - (1 - s)i/"){l + a(l - - 

where n(t, s) is a product copula of independent random variables. This means 
that if the joint distribution of {X,Y) is FGM given in (|23p . then the extreme 
order statistics Xn.n and are asymptotically independent. So are the order 
statistics Xi-n and Yim. Despite of this fact, the Baker's type BB distribution 
constructed on the base of the copula (|24p and (P5|) has enough large correlation. 
Indeed, consider Baker's type BB distribution © with the "kernel" copula (|24p 
and with the "kernel" copula ([25]) : 



1 " 

n ^ — ' 

k=l 
1 ^ 

1 " 

-ECi,„.„(F|"(x),F^^"(y)) 



/c=l 



1 " 

n ^-^ L 



- 1 



(1 - + (1 - - 1 



+ (1 - (1 - - (1 - 

x{l + a(l - - "} . 

For example, the correlation coefficient of G^'"* (a;, y) for Uniform{0, 1) marginals 
is = 0.4985 for n = 2. 



3.2 The Baker's type BB distributions with "kernel" cop- 
ula of bivariate Gumbel's extreme order statistics 



Let 



Ciu,v)^— . (26) 

1 + u — uv 



is the copula of the Gumbel's bivariate logistic distribution function 

Hx,Yix, y) = {l + e-" + e-y)-\ x>0,y>0 
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with standard logistic marginal distributions Hx{x) = {1+e ^) ^ and Hyiy) = 
(1 + e-f)-\ X >0,y >0. (Gumbel (1961), see also Nelsen (2005), page 28). 
Using 



Fxix)+FYiy)-Fxix)FY{y)' 



one obtains from (pUj) 

Cn:n (^; •^) 

and from (I^TI) 



Cl:„(t,s) = t + S-l + 

(l-i)i/" + (l-s)i/"-l 

(l-(l-t)l/«)(l-(l-5)l/n) 



(27) 



(28) 



2 - (1 - f)i/" - (1 - s)i/" - (1 - (1 - t)i/")(l - (1 - s)i/") 



It is seen that 

Jim C„:„(i, S) = Jim (,V»+,l/n*l^l/n,l/»)n = ^ 

4 Joint distribution of bivariate order statistics 
for Frechet upper bound copula 

Let {Xi, Yi), {X2, Y2), {Xn, F„) be a bivariate sample with joint distribution 
function F{x,y) = C {Fx (x) , Fy (y)) ■ In this section our aim is firstly, to 
investigate the joint distribution function of bivariate order statistics {Xr-.n, Yr-.n) 
for a copula with maximal correlation, i.e. the Frecheet upper bound. We are 
interested then in distribution function 

-. n n n b 
r— 1 i—r j—r k—a 

in the case where the marginal distributions are uniform and C(w, v) =^ min(u, v). 
Recall that the coefficients c(n, k; and pii,pi2,?'2i,P22 are given in ([H]),®. 
Secondly, we consider a distribution introduced in ([9]) with uniform marginals, 
i.e. 

n 

G^+i^^y) = -Y.C{F^x''{x),F^--{y)) 



n 



fc=i 
1 " 

-T.GiFu.:A^),Fv_{y)), 
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where C{t,s) = mm{t,s) and Fu^.^{x) = ^ - xy-' and Fv,,^{y) = 

i—r 

i—r 

Assume that marginal distributions are Uniform{0, 1), i.e. Fx{x) — x,0 < 
X < 1, Fyiy) = y, < ?/ < 1. Denote by {Ui,Vi),i = l,2,...,n the random 
sample from the bivariate distribution C{u,v), < u < 1,0 < v < 1 and 
Ui:n < U2:n < ' ' ' < Un:n, Vi-.n < V2:n < ' " ' < K:« be Corresponding order 
statistics. Then from (jl]) for r = s one has 

P{Ur:n < U, Vr:n < «} 
n n b 

i—r j— r k—a 

x{v- Ciu, v)y-''iC{u, 0<w<l,0<v<l. (29) 

and 

n n n b 
r— 1 i—r j—r k—a 

[v - C{u,v)y-^{C{u,v))'^-"'-^+^ ,0 < u < 1,0 < < 1. (30) 
Lemma 1. The joint distribution of Ur-.n and Vr-.n can be represented as 

P{Ur:n < U.Vr.n < v} 

= E f") C\u,v)[v-CM+C{u, v)] " 



+ C-'' (m, i;) [li - C(u, v) + C(?/, v)] " ' 

j=r 
n n 

+ EE E c{n,k;i,j)C{u,v)Hu-C{u,v)y-'' 

i—r j—r k^i^j 

x{v- Ciu, v)y-''iCiu, v))"-'-^+'=. (31) 
Proof. Separating terms for summation in (j29p for k = i=j,k = i=/=j and 
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k = j i we have 

P{Ur:n < U, K:n < v} 

n 

= c{n, i, i, i)C\u, v){C{u, v)Y 



n n 



n n 



+ ^ ^ c(n, i, i, j)C^ {u, v){u - C{u, v)Y-^ [v - C(u, v)y-^ {C{u, 

n n 

+ E E E ^("' ^; *'-?')^("' - ^("' «))'^'(« - ^("' ^))'"'(^(«' 



n 
i 

n n 



= E(")^^("'^)(^("'"))' 



7,1 



- - j)! 
n! 

j\{i-j)Hn - i)! 



C\u, v){v - C(m, v)y-'{C{u, v))"-^' 

c^' {u, v){u - ciu, v)y-^ {ciu, v)y'-' 



+EE 

i—r 
n n 

+EE 

n n 

+ EE E cir>;k;i,j)C{u,vy{u-C{u,vy 

i—r j—r k^i^j 

x{v- Ciu, vy^''' iC{u, 



i — k 



(32) 



Consider the second summation in p2p and changing index in the inner sum as 

j ^ i — k {j — i + k, j = i + l k ~ 1, j — n k = n — i) we have 



EE 

i—r j—i-}-! 
n 

J2C\u,v) 



j)! 



c\u, v)iv - c{u, vy^-'ic{u, vy 



E 



{n — iy. 



i\{n — iy. ^-^ kl{n — i — k) 



■{v-ciu,vy\ciu,vy 



I— i — k 



E 

.fc=0 



(n-iy 



k\{n — i — ky. 



{v-ciu,vy''{c{u,vy 



71— i — k 



(C(u,«))"- 



J2(^^c\^,y)[y-ciu,v) + c{u,v)]'' ' -j2{^^c\u,vyc{u,vy^^-\ (33) 
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Analogously, the third term in (15^ can be written as 

n n I 

E E _ - .y C'Miu - C{u,v)y-^iCiu,v)r~^ 

j=r i=j + l -^'^ •'''^ 

= izi'^^^'Mb-CM+CMY-' (34) 



j=r 



Taking into account ^ and ^ in ([22]) we have (I5T|) . 
Theorem 2. If C(m,w) = min(u,w), then 

P{Ur:n < U, Vr:n < v} 
n 

j2{^i)u^[i~ur-' u<v 

i—r 
n 

i—r 

/ P{Ur:n 1^ u) if U < V 
{ P{Vr:n < v} if U>V 

Proof. Consider (PT|) . Let C{u,v) = niin(u,w). Then it is clear that 

(u — C{u,v)y~'^{v — C{u,v)y~^ = for those k satisfying k ^ i and k ^ j, 
because u — C{u, v) = u — u — Qii u < v and v — C{u, u) = if u > v. Therefore 
the last term of (1311) vanishes and we have from (|3T]) 



(35) 



P{Ur:n < U, Vr:n < v} 

f^(^)uMl-<-' if u<v 

i—r 
n 

5](^)i;ni-<"^ if u>v 

i—r 

J P{Ur:7i < u} if U < V 
~ \ P{Vr:„ <V} if U> V 

Consider now, the Baker's Type I BB copula obtained from ([35 

1 " 

{U, v) = -y P{Ur:n < U, Vr:n < v} 
r=l 

n n 

^ EE u<v 

r—l i—r 
n 71 

^EE(:)-Mi-<"' If u>v 



r— 1 i—r 
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if u <v 



= W{u, v) = niin(u, v) 



V if u > V 

is nothing but the Prechet upper bound itself. In other words the Baker's 
BB distribution with uniform marginals and underlying joint distribution being 
Frechet upper bound generates the Frechet upper bound. 

Now consider the Baker's Type II BB distribution with the "kernel" copula 
being Frechet upper bound, i.e. 

1 " 
n 

r=l 

n 

where C{t,s) = min(i,s) and Fu,.^{x) = ^ (")a;*(l - x)"^* and Fy^^Jy) = 

i=r 

n 



It follows that 



1 

r=l 

1 " 

-Y,^HFur.Ax),Fv^.,Jy)) 

n 
n 

^E"=iE(")2^'(i-2/)""*) ^>y 

i=r 

= m.m{x,y). 

Therefore, if in the constructions K^\x, y) and G^]^\x, y) with Uniform{0, 1) 

marginals one uses the Frechet upper bound copula W{f,, s), then obtained cop- 
ula is again the Frecheet upper bound copula. From the Example 1 one ob- 
serves that if in the same construction one uses FGM copula then different 
distributions with high correlation can be obtained. For different constructions 
satisfying conditions of the Example 1, we have the following set of inequalities 
for correlation coefficient 

Pm < Pg^^) ^ /'ifi") ^ Ph^^> - Ph'-^^ - Pk'-^'^ - Pg'-:^^ - Pw' 

where pj^f is the correlation coefficient of the Frecheet lower bound M{t, s) = 
max.{t + s — 1, 0). 
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